#include "mex.h"

static void mycalc(double p[], double n[], int r, int c)
{
        int i;
        for (i=0; i<r*c; i++)
            p[i]=2;
        
        p[0]=1;
}

void mexFunction(int nlhs,mxArray *plhs[], int nrhs, const mxArray *prhs[])
{
        double *p, *n;
        int r,c;
        r=mxGetM(prhs[0]); c=mxGetN(prhs[0]);
        plhs[0]=mxCreateDoubleMatrix(r,c,mxREAL);
        p=mxGetPr(plhs[0]);
        n=mxGetPr(prhs[0]);

        mycalc(p,n,r,c);
}
